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ABSTRACT 



Context. It is still debated whether star formation process depends on environment. In particular it is yet unclear whether star formation 
in the outer Galaxy, where the environmental conditions are, theoretically, less conducive, occurs in the same way as in the inner 
Galaxy. 

Aims. We investigate the population of NGC1893, a young cluster (~ 3 - 4 Myr) in the outer part of the Galaxy (Re > 1 1 Kpc), to 
explore the effects of environmental conditions on star forming regions. 

Methods. We present infrared observations acquired using the IRAC camera onboard the Spitzer Space Telescope and analyze the 
color-color diagrams to establish the membership of stars with excesses. We also merge this information with that obtained from 
Chandra ACIS-I observations, to identify the Class III population. 

Results. We find that the cluster is very rich, with 242 PMS Classical T-Tauri stars and 7 Class 0/1 stars. We identify 110 Class III 
candidate cluster members in the ACIS-I field of view. We estimate a disk fraction for NGC1893 of about 67%, similar to fractions 
calculated for nearby star forming regions of the same age. 

Conclusions. Although environmental conditions are unfavorable, star formation can clearly be very successful in the outer Galaxy, 
allowing creation of a very rich cluster like NGC1893. 

Key words. Stars: pre-main sequence, formation, circumstellar matter. Open cluster and association: individual: NGC1893 



1. Introduction 

A key issue of the star formation process is its dependence on en- 
vironmental conditions. Theoretically, temporal and spatial vari- 
ations are expected since, for example, a fundamental parame- 
ter of star formation, the Jeans mass, depen ds on gas temper- 
ature, chemical abundances, and density (cf. lElmegreenl I2002L 
2004). Moreover the variation in the initial mass function at low 
stellar masses across different environments suggests that cloud 
fragmentation based on the Jeans formulation, i.e., on gravity, 
is not the only mechanism at wor k, but that other processes 
(e.g. turbulence) are also im portant dPadoan & Nordlundl 1999; 
iMac Low & Klessenlr2004l) . 

We observe stars forming under very different conditions and 
one of the key questions is whether star formation proceeds in 
the same way under the very different environmental conditions 
of the inner and outer region of the Galaxy. Molecular clouds 
and very young stellar associations are also present a t large dis- 
tances from the Galactic Center (ISnell et all Hp02), implying 
that star formation still occurs in these regions. Conditions in 
the outer Galaxy should be much less conducive to star forma- 
tion: the average surface and volume densities of atomic and 
molecular hydrogen are much smaller than in the sola r neigh- 
borhood or in the inner Galaxy ( Wouterloot et al l 1 19901) . the in- 
terstellar radiation field is weaker dMathis et all Il983l) . promi- 
nent spiral arms are lacking and there are fewer supernovae to 



* Tables [T] and [4] are only available in electronic form at the 
CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat7J/A-l-A/ 



act as external trig gers of star format i on. M etal content is, on 
average, smaller (Wilson & Matteucci. [19921) . decreasing radia- 
tive losses and therefore increasing cloud temperatures and con- 
sequently pressure support . Moreover, the pr essure of the inter- 
cloud medium is smaller (Elmegreen, 1989). Given these rel- 
evant differences, quantitative observations of star forming re- 
gions in the outer Galaxy may help to establish the role of phys- 
ical conditions in the star formation process. 

In order to investigate the star formation processes in the 
outer Galaxy, studied the young cluster NGC1893, whose galac- 
tocentric distance is > 1 1 Kpc. Lying in the Aur OB2 association 
toward the Galactic anti-center, NGC1893 is associated with the 
HII region IC 410. It contains a group of early-typ e stars with 
some m olec ular clouds but o nly m oderate extinction; Tapi a et all 
(11991b and IVallenari etaO ([1999) estimated an age of 4 Myr 
from stars already on the Main Sequence, but the presence of 
two Q5 V stars indicates that the regio n should be younger than 
~3 Myr; iMarco & Negueruelal (|2002|) found PMS B-type stars 
consistent with age < 1 Myr. These results suggest that star for- 
mation is still ongoing and hence there may be a range of ages 
for stars in the cluster. 

Several authors performed UBV photometry of the clus - 
ter (ICuffevL fT97l iMoffat & VoetL fl974t iMassev et all fl995l) . 
IVallenari et all (1 19991) obtained NIR photometry of the cluster, 
concluding that NGC1893 ma y have a rich pre-main sequence 
population. iTapia et all (119911) determined an Ay = 1.68 and a 
dereddened distance modulus of (V - My)o = 13. 18 + 0. 11, cor- 
responding to a heliocentric distance of 4.3 Kpc at / = 173.6°; 
other authors calculated the distance modulus corresponding to 
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a distance in a range between 3.2 and 6 Kpc ( Massev et al 



Sharma et al 



199 1 iMarco etalL l200lt iDaflon & CunhaL 120041 
2007). With the solar galactocentric distance taken to be 8.5 
Kpc, NGC1893 is therefore located at a distance of 11.7 - 14.5 
Kpc from the Galactic Center. NGC1893 has a s light! 
metal content, as reported by Rolleston et ak| d20Qi 
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analysis of spectra of a few B stars. IDaflon & Cunhal (12004) con- 
firmed that the abundances of C, N, and O in NGC1893 were two 
to three times lower than in the Sun. 

We present the first results of the joint Chandra-Spitzer large 
project The Initial Mass Function in the Outer Galaxy: the star 
forming region NGC1893 aimed at studying star formation pro- 
cesses at high galactocentric distances using multi-wavelength 
data, including infrared, X-ray, and optical observations. By 
means of the infrared observations, we study the properties of 
members with disks and/or envelopes, while the analysis of X- 
ray emission is crucial to establishing the membership of stars 
without infrared excesses and is also the elective way to inves- 
tigate the coronal properties. By means of the optical observa- 
tions, we finally estimate the masses of the stars and determine 
the distance and the age of the cluster. The project, led by G. 
Micela, includes a ~ 3 1 .2' x 26.0' IRAC/Spitzer map of the clus- 
ter and 5 Chandra ACIS observations with a total exposure time 
of ~ 440 Ks. 

In this first work, we conduct a census of the pre-main sequence 
candidate members of NGC1893. The IRAC observations of 
NGC1893 allow us to determine the members of the cluster with 
infrared excesses: Class 0/1 and Class II stars, while the 5 ACIS 
observations are used to select the Class III objects that have 
already lost their disks and have no prominent infrared excess. 

In Sect. [2] we describe the IRAC observations and the data 
reduction. In Sect.[3]we describe the X-ray observations, in Sect. 
|4] we identify the members of the cluster, while in Sect. [5] we 
summarize our results. 

2. Infrared Observations and Data Reduction 

2.1. IRAC Observations 

NGC1893 was ob served with the Infrared Array Camera (IRAC; 
Fazi o et ail |2004|) onboard the Spitzer Space Telescope on 2006 
March 25. This instrument has four detector arrays, with each 
dedicated to a different wavelength band, centered at 3.6, 4.5, 5.8 
and 8.0 pm, each of which has a field of view of ~ 5.2' x 5.2'. 
All four bands were observed simultaneously, using the High 
Dynamic Range (HDR) mode with a long integration frame of 
26.8 s and a short one of 1.0 s. We covered the region of the 
cluster by a 5 rows x 5 columns map with five dithers at each 
map position. The total observed area was ~ 31.2' X 26.0'. 

The data reduction started from the IRAC Basic Calibrated 
Data (BCD), provided by the Spitzer Science Center (SSC) 
pipeline reduction (version S13.2.0). The SSC process includes 
flat field corrections, dark subtraction, linearity, and flux calibra- 
tion. The pipeline also masks instrumental signatures, bad pix- 
els, and saturated pixels due to bright sources. 

We further processed the BCDs by means of the IRAC arti- 
fact correction code, in order to mask bad pixels not identified 
by the SSC pipeline reduction, and correct muxbleed and column 
ull-down defects occurring close to moderately bright sources 



Using the SSC MOPEX code dMakovoz & Marleaul 12005b . 
we combined the BCDs into a single mosaic image for each 

1 All the documentation for this code is available at http:// spi- 
der, ipac. caltech. edu/staff/carey/irac ^artifacts 



channel and exposure time, by averaging together the overlap- 
ping frames; in this way, we obtained four images for the long 
integration frame time and four for the short one, with a S/N 
factor higher than for the individual BCDs. 

Figure Q] shows the mosaics obtained for the 4.5 and 8.0 pm 
channels. The figure illustrates that the two short-wavelength 
channels are far more sensitive for detecting stars than the two 
longer wavelength channels both because the detectors for the 
3.6 and 4.5 pm channels are far more sensitive and because nor- 
mal stellar photospheres are brighter at 4.5 than 8.0 microns. The 
nebula associated with NGC1893, however, is far more evident 
at 8.0 pm because there are prominent PAH emission features in 
that band. 



2.2. Infrared Photometry and Catalogue 

We an alyzed each m osaic image with the DAOPHOT II proce- 
dures dStetsonl 1 1987b . to obtain a list of point sources, their po- 
sitions, and their magnitudes. Since DAOPHOT expects images 
calibrated in electrons (or electrons/second), we converted the 
surface flux images back into counts, using the flux conversion 
values found in the IRAC Data Handbook We found star-like 
objects using a detection threshold of 5<x for 3.6 and 4.5 pm 
channels and 6cr for 5.8 and 8.0 pm channels, where cr is the 
standard deviation of sky counts 0. 

We derived aperture photometry with an aperture radius of 3 
pixels (3.6"), considering an annulus with inner radius 12" and 
outer radius 24" for the background extraction, and we applied 
the relevant aperture corrections. For the zero point fluxes and 
the aperture correcti ons, we used values found in the IRAC Data 
Handbook (s ee also iReach et all 120051) . Using the DAOPHOT 
II/ALLSTAR dStetsonl 1 1987b and ALLFRAME (IStetsonL [1994) 
procedures, we also computed the PSF photometry, calibrating 
the PSF magnitudes by means of a comparison with the aperture 
photometry. Since, for the bright stars, the uncertainties related 
to the PSF magnitudes are larger than the aperture ones, and, 
on the other hand, the PSF errors are small for faint stars, we 
decided to use, for each star the magnitude computed with the 
procedure providing the smallest uncertainty. 

In order to eliminate false sources, we merged the four 
IRAC source lists and the 2MASS K band catalogue, by cross- 
identifying sources, using a radius of 0.6". We accepted a star in 
the final catalogue if we were able to identify it in at least two 
different wavelengths. Moreover we considered only the magni- 
tudes that have an uncertainty of less than 0.1. In Table Q] we 
provide our catalogue, which includes the sequential numbers, 
the celestial coordinates, the magnitudes in K, and in 3.6, 4.5, 
5.8, and 8.0 pm and the DAOPHOT sharp parameter, that mea- 
sures how well the brightness distribution could be interpreted 
as a point-like function (we use this parameter in the following 
to identify extended sources). 

In Table 12 we report the total number of detected objects in 
the catalogue, specifying the number of sources in each IRAC 
channel and in the 2MASS catalogue and the detection limit for 
each band. At 5.8 pm and 8.0 pm, where the instrument sensitiv- 
ity is lower than in the other two channels and the background is 
higher, we find fewer sources. Although the number of sources 
common to the four channels is strongly limited by the sensitiv- 
ity of these two channels, we found 1028 objects, 25% of which 
were not previously detected with 2MASS. 



2 http:// ssc. spitzer. caltech. edu/irac/dh/iracdatahandbook3 . 0.pdf 

3 The detection threshold is higher in 5.8 and 8.0 [im channels to 
avoid spurious sources corresponding to the nebula peaks 



Table 1. Data Catalogue. ID is the source identification number, RA(J2000) and DEC(J2000) are the celestial coordinates, K is the 2MASS K magnitude and <r K is its 
photometric uncertainty, [3.6], [4.5], [5.8], [8.0] are the magnitudes in the four IRAC channels and crp.6] , C[4.5], C[5.8]> C[8.0] are their photometric uncertainties. The complete 
table is available in electronic format at the CDS. 



ID 


RA(J200Q) 


D£C(72000) 


K 


<r K 


[3.6] 


Cr [3.6] 


[4.5] 


(T [4.5] 


[5.8] 


cr [5.8] 


[8.0] 


cr [8.0] 


sharp 


26 


5 


: 21 


51.59 


33 : 


30 : 


43.03 


10.62 


0.02 






10.492 


0.003 


11.006 


0.095 


10.603 


0.071 


-0.486 


91 


5 


: 21 


52.94 


33 : 


21 


49.24 


10.60 


0.02 






10.506 


0.002 


10.556 


0.048 


10.389 


0.043 


0.065 


121 


5 


: 21 


53.43 


33 


: 27 


: 9.55 


12.31 


0.02 


11.367 


0.003 


10.865 


0.004 


10.230 


0.004 


9.700 


0.008 


-0.001 


142 


5 


: 21 


53.74 


33 : 


23 


26.93 


13.01 


0.03 


12.777 


0.013 


12.748 


0.011 


12.041 


0.081 






0.234 


151 


5 


: 21 


53.87 


33 


: 26 


: 8.63 


9.64 


0.02 


9.441 


0.046 


9.397 


0.081 


9.413 


0.004 


9.378 


0.016 


-0.027 


159 


5 


: 21 


53.96 


33 : 


34 : 


55.75 


12.60 


0.02 


12.391 


0.006 


12.424 


0.008 


12.421 


0.039 


12.592 


0.098 


-0.082 


166 


5 


: 21 


54.02 


33 : 


31 


42.45 


13.58 


0.03 


13.685 


0.031 


13.629 


0.031 


13.609 


0.041 






0.015 


176 


5 


: 21 


54.15 


33 : 


34 : 


24.26 






13.811 


0.069 






11.413 


0.077 


9.602 


0.090 


0.664 


183 


5 


: 21 


54.26 


33 


: 35 


: 6.90 


14.40 


0.06 


14.691 


0.042 


14.298 


0.034 


14.069 


0.091 






0.185 


189 


5 


: 21 


54.30 


33 : 


27 : 


42.04 


13.97 


0.05 


13.699 


0.028 


13.742 


0.036 


13.513 


0.058 






0.160 


196 


5 


: 21 


54.38 


33 : 


27 : 


34.23 






15.924 


0.044 


15.082 


0.057 


14.164 


0.082 






0.067 


197 


5 


: 21 


54.43 


33 


: 32 


: 0.52 


12.36 


0.02 


12.326 


0.007 


12.290 


0.008 


12.046 


0.041 


11.280 


0.079 


0.156 


236 


5 


: 21 


54.85 


33 


: 27 


: 5.93 


14.71 


0.08 


14.499 


0.029 


14.523 


0.033 


14.069 


0.085 






0.138 


238 


5 


: 21 


54.89 


33 


: 33 


: 9.28 


12.60 


0.02 


12.413 


0.005 


12.530 


0.006 


12.452 


0.028 


12.362 


0.094 


0.022 


246 


5 


: 21 


55.00 


33 : 


30 : 


50.75 


11.87 


0.02 


11.682 


0.011 


11.827 


0.013 


11.860 


0.063 






-0.109 


297 


5 


: 21 


55.66 


33 : 


29 : 


23.51 






13.684 


0.083 






11.179 


0.093 


9.277 


0.092 


0.576 


301 


5 


: 21 


55.75 


33 : 


22 : 


42.11 


12.24 


0.02 


12.206 


0.009 


12.068 


0.008 


12.189 


0.025 


12.106 


0.043 


0.115 


302 


5 


: 21 


55.76 


33 : 


33 : 


52.78 


12.22 


0.02 


11.903 


0.005 


11.985 


0.007 


11.950 


0.026 


11.817 


0.097 


0.060 


303 


5 


: 21 


55.79 


33 : 


26 : 


34.29 






14.866 


0.031 


15.203 


0.081 


12.423 


0.044 


10.591 


0.035 


0.776 


305 


5 


: 21 


55.82 


33 


: 35 


: 8.31 


13.71 


0.04 


13.699 


0.011 


13.773 


0.012 


13.635 


0.098 


13.159 


0.086 


0.141 


325 


5 


: 21 


56.05 


33 : 


29 : 


27.15 






13.597 


0.081 






11.138 


0.097 


9.245 


0.094 


0.592 
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Fig. 1. Images of the ~ 31.2' x 26.0' mosaic for IRAC 4.5 (left) and 8.0 pm (right) channels. 



Table 2. Number of infrared detected sources in the NGC1893 region and detection limits 



list Number of sources Detection limit (mag) 



All catalog * 15277 

3.6 yum 12176 19.03 

4.5 urn 10316 18.34 

5.8 yum 2932 15.88 

8.0 yum 1705 14.92 

K2MASS 2203 15.01 
4 IRAC bands 1028 
4 IRAC bands and K 2MASS 77 1 



* Sources with at least two magnitudes and for which at least one magnitude has photometric error < 0.1. 



3. X-ray Observations and Data Reduction 

The X-ray observations of NGC1893 combine four nearly con- 
secutive exposures of the cluster taken in 2006 November and 
a fifth exposure taken in 2007 January (see Table [3] of details), 
with a total exposure time of ~ 440 Ks. The detailed analysis of 
these X-ray data will be the subject of a future paper. Here we 
use the X-ray observations to identify probable members among 
the IRAC sources without any infrared excess in order to es- 
timate the number of Weak-lined T Tauri stars. Therefore, we 
summarize only the main steps of the analysis. 

The observa tions were obtained with the ACIS camera o n- 
board Chandra (Weisskopf et al., 2002t iGarmire et all 120031) in 
the FAINT mode; we considered only results from ACIS-I, the 
four front-side illuminated CCDs covering about 17' x 17' on the 
sky. ACIS-I is characterized by high angular resolution, with an 
on-axis PSF FWHM of ~ 0.5", which makes it possible to re- 
solve the central region of the cluster, despite the large distance 
ofNGC1893. 

We detect ed sources applying the PWDetect code 

dDamiani et aU 1 19971) to the merged event file of the five 



observations, filtering photons with energy between 0.5 and 8.0 
KeV. The significance threshold was set to 4.6 cr, corresponding 
to 10 spurious sources detected in source-free fields with the 
background level of our observations. PWDetect reported 1025 
sources. Upon careful inspection, we removed four entries 
corresponding to sources that were detected twice, leaving a 
total of 1021 distinct sources. 



In order to identify young stars in which the primordial disk 
had already dissipated, we combined the Spitzer and Chandra 
data. Since the field of view of Chandra ACIS-I is limited to 
17' x 17', we were unable to find X-ray counterparts for all ob- 
jects in the Spitzer field of view. For comparison, we overplot 
the ACIS field of view on the IRAC image at 4.5pm (Fig.©. 



Due to the degrading of the ACIS PSF at high off axis, we 
used three different matching radii (0.7", 1.5", 2.5"), depending 
on the source off axis angle. The results of the match between 
the IR and the X-ray catalogues are given in Sect. 14.1.31 
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Table 3. X-ray Chandra ACIS observations 



ObsID 


Start Time 


Exposure 


R.A. (J2000.0) 


Deck (J2000.0) 


Roll Angle 






(ks) 


(deg) 


(deg) 


(deg) 


8462 


2006 Nov 07 13:33:15 


42.61 


80.7088633 


33.4726607 


107.50840 


6406 


2006 Nov 09 12:51:18 


115.7 


80.7088706 


33.4726628 


107.50839 


6407 


2006 Nov 15 05:31:19 


126.2 


80.7088692 


33.4726604 


107.49471 


8476 


2006 Nov 17 10:55:16 


53.26 


80.7088737 


33.4726629 


107.49472 


6408 


2007 Jan 23 00:12:30 


102.8 


80.7054833 


33.4640840 


262.11412 




Fig. 2. ACIS field of view (17' x 17') on the IRAC image at 4.5 fjm. 



4. Young Stellar Objects in NGC1893 

4.1. IRAC four band source characterization 

In this section we analyze the subsample of sources detected in 
all four IRAC channels to find the likely members with excesses. 
We first reject the contaminants (extragalactic sources or nebular 
emission) and then classify the young stars with excesses. 

4.1 .1 . Estimate of Contaminants 

Young stars with disks are characterized by infrared excesses, 
therefore the study of their Spectral Emission Distribution (SED) 
(either via their SED slope or their position in color-color dia- 
grams) is the most powerful method for distinguishing the clus- 
ter members with excesses from the old stars in the same field 
of view. Nevertheless, there are at least two classes of con- 
taminants: the extragalactic sources, mainly broad-line AGN 
(BLAGN) and star forming galaxies, and emission from the neb- 
ula that, for very embedded stars could contaminate the photom- 
etry. 



To avoid c ontamination f rom e xt ragalactic sourc es , sev- 
eral authors dHarvev et ail I2006L |j0rgensen et al.L [2006; 
IHernandez et al.L I2007I) in their study of nearby star forming 
regions, suggested a cut in magnitude at 3.6 /im or at 4.5 jum. 
However, this method alone is less effective in the case of 
NGC1893 because it is much further from the Sun than the star 
forming regions considered in those studies, and its PMS stars 
are therefore faint. 

Following IStern et al.l d2005l) and iGutermuth eTail d2007l) . 
we decided to identify star-forming galaxies by means of their 
typical Polycyclic Aromatic Hydrocarbons (PAH) emission. 
Using the [4.5]-[5.8] vs. [5.8]-[8.0] and the [3.6]-[5.8] vs. [4.5]- 
[8.0] color-color diagra ms in Fig. [3] and the equations given in 
IGutermuth et all (120071) . it is possible to distinguish the regions 
delimited by the solid lines in Fig. [3] where objects whose SED 
are dominated by PAH emission (indicated by empty squares) 
are typically found. 

Since only a small number of young stars were shown to 
have significant PAH emission in other star forming regions 
(lHartmann et al.L 120051: IHernandez et al.L l2007h . we rejected the 
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Fig. 3. IRAC color-color diagrams illustra t ing the selection of star-forming galaxies, characterized by a strong PAH emission. Using 
the method described in Gutermu th et al.l ([2007), we rejected the sources lying in the region delimited by the solid line in at least 
one of the two color-color diagrams. As explained in the legend, we designate the PAH-rich galaxies with little squares. 



32 sources that satisfied at least one of the conditions described 
ld2007l) . 



in Guterm uth et al] 

IStern et ail d2005l) found broad-line AGNs in the following 
region of the [3.6]-[4.5] vs. [5.8]-[8.0] diagram: 



[5.8] - [8.0] > 0.6 
[3.6] - [4.5] > 0.2 • 
[3.6] - [4.5] > 2.5 • 



([5.8] -[8.0]) + 0.18 
([5.8] - [8.0]) - 3.5 



(1) 



Due to the distance of NGC1893, a simple cut in magnitude 
would be ineffective to filter AGNs alon e, and, on the other 
hand, the zone found by [Stern et al] (120051) is similar to that oc- 
cupied by protostars; we therefore decided to combine the two 
conditions, selecting AGNs as the sources that follow the con- 
ditions CD and whose 3 .6 /mi magnitude is greater than 14.5 
dHernandez et all 120071) . Using this criterion, we identified 56 
candidate AGN. 

As sources contaminated by the nebular emission, we se- 
lected 23 objects with sharp dStetsonl [l987h greater than 0.37, 
corresponding to 3cr of the sharp distribution , i.e. objects with 
a brightness distribution inconsistent with a point-like source. 
These sources have very red [5.8]-[8.0] colors, and most have a 
blue [3.6]-[4.5] color, due to the fact that the nebula is brighter 
at 4.5 /im. It is possible that some of these objects are embedded 
members of the cluster and we do not classify them because their 
photometry could be contaminated by the nebular emission. 

We illustrate the contaminants we identified by means of 
these criteria in Fig. |4] Squares represent the candidate star- 
forming galaxies, found by means of diagrams in Fig. [3] crosses 
(above the solid line) indicate candidate AGNs, and triangles 
represent sources contaminated by the emission of the nebula. 

4.1.2. Members with infrared excesses 

There are two standard ways to classify young stellar o bjects, 
both based on the analysis of their SED. The first method (lLadal 



I1987I) uses the shape of the SED, which in IRAC bands can be 
approximated to a po wer law. From a co mparison between mod- 
els and observed SED. lLada et al.l(l2006l) found that, in the IRAC 
bands, stars with primordial disks have a slope a > —1.8; stars 
with purely photospheric emission should have a < -2.56; stars 
with slopes intermediate between these two values are plausi- 
bly "transitional objects" that still have disks but possibly with 
inner holes or with optically t hin emission regions. Another stan- 
dard method in the literature (lAllen et all 12004; Me geath et all 
2004) consists of the conversion of SED models into position in 
the [3.6]-[4.5] versus [5.8]-[8.0] diagram. 

Here we used a method similar to that o f lAllen et all (12004) 
and iMegeath et al.l d2004l). with an extensive grid of models of 
young stars (Robitaille et al., 2006), to define regions in color- 
color diagrams rela tive to different P MS cl asses. 

Following the Robitai He et al] (j2006) notation for which 
M e „v is the accretion rate from the envelope, M+ is the mass of 
the central star, and M^isk is the mass of the disk, it is possible to 
define class 0/1 stars as those with: 



M mv /M* > 10- 6 yr" 



and class II stars as: 
yr 



M CTV /M* < 10- 6 



M^t/M* > 10 



-6 



(2) 



(3) 



Converting the SEDs of th ese groups of models in to IRAC colors 
(see figures 17 and 18 in iRobitaille etafl 2006). we note that 
the class II objects occupy the region of the color-color diagram 
defined by [3.6] - [4.5] < 0.8 and [5.8] - [8.0] > 0.40. 



4 We analyzed only models with inner disk radius less than 100 times 
the dust sublimation radius. Stars with internal disk radius greater than 
this limits have [3.6]-[4.5] and [5.8]-[8.0] colors close to zero, and the 
first color different from zero is [8.0]-[24.0], where we do not have ob- 
servations. 
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Fig. 4. IRAC color-color diagram highlighting the contaminants 
to the young star sample. As explained in the legend, squares 
indicate the star-forming galaxies selected by means of the di- 
agrams in Fig. [3] crosses indicate the BLAGN, selected as the 
weak sources above the solid line, while triangles indicate the 
sources contaminated by the emission of the nebula: most of 
them are characterized by a strong red color in [5.8]-[8.0] and 
a blue color in [3.6]-[4.5]. 



Analyzing the colors relative to class 0/1 protostars, the 
models occupy a large region of the color-color diagram be- 
cause of temperature and inclination effects. Although the locus 
[3.6] - [4.5] > 0.8 is typical of class 0/1 objects, the region occu- 
pied by class II stars can also be occupied by protostar models. 
Therefore, although many class 0/1 stars can be uniquely iden- 
tified, the evolutionary classifications of the remaining class 0/1 
as well as most class II objects cannot be derived from the IRAC 
colors alone. However, considering the different duration of the 
two stages, in a cluster of 3 Myrs it is 10-100 times less prob- 
able to find a class 0/1 than a class II star. Therefore, we des- 
ignate as class 0/1 only the 7 stars whose colors place them in 
the region uniquely corresponding to protostars, while we desig- 
nate the 242 stars lying in the color-color diagram region com- 
mon to the two evolutionary stages as class II objects. Using 
this method, we are able to identify 249 members of the cluster. 
In Fig. [5] we plot the IRAC color-color diagram for our stars, 
using the previous classification. As explained in the legend, tri- 
angles represent the class 0/1 objects, while the crosses indicate 
the class II stars. 



4.1 .3. Weak lined T Tauri 

Although the star formation process is often considered a low 
energy process, high energy processes occur as early as the pro- 
tostellar phase and pre-main sequence (PMS) stars are strong X- 



Fig.5. [5.8]-[8.0]vs.[3.6]-[4.5] IRAC color-color diagram. 
Triangles represent the objects whose SED is dominated by the 
envelope emission, while crosses refer to objects with SED dom- 
inated by the disk emission. The arrow refers to the reddening 
vector for Ay = 10. To better classify the stars at the border 
between these two classes, we took into account the error asso- 
ciated with the colors. 



ray sources with luminosity up to 10 times the X-ray luminosity 
of the present-day Sun. Among the sources detected by IRAC, 
we utilize this property to identify young stars that have already 
dissipated their disks and would otherwise be indistinguishable 
from field stars (Class III or Weak lined T Tauri). 

Since we are able to compare only the objects in the same 
field of view, we consider for this analysis the 12513 infrared 
sources in the 17'x 17' AClS-l/Chandra FOV. Applying the anal- 
ysis described in Sect. 14.1. Ti to the 692 sources with emission in 
all four IRAC channels, we labeled 77 sources as possible con- 
taminants. As discussed in Sect. 14.1.21 by means of the IRAC 
color-color diagram we identified 219 stars with disks and 6 ob- 
jects with envelope emission. 

In our ACIS observations, we identify 1021 X-ray sources. We 
expect that this sample is strongly contaminated by extragalactic 
sources because NGC1893 is located towards the anticenter di- 
rection. Using the galactic H column density Nh ~ 6 x 10 21 crrT 2 
0, the flux limit of our X-ray obs ervations and the results from 
the Chandra North Deep Field (Bran dtet all 1200 It) , we predict 
300-550 extragalactic contaminants in our X-ray catalogue. 

By matching the X-ray catalogue with all the 12513 infrared 
sources in our catalogue of sources in the ACIS field of view, we 
find 720 matches (~ 70% of the X-ray sources). Considering the 
692 infrared sources detected in all the four IRAC energy bands, 
we find 248 matches, 20 of which correspond to extragalactic 
contaminants and one to a source contaminated by the emission 



5 calculated by means of the Nh tools of HEASOFT (see 
http://heasarc.gsfc.nasa.gov/cgi-bin/Tools/w3nh/w3nh.pl ) 
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Fig. 6. IRAC color-color diagram for stars in the Chandra ACIS- 
I field of view. Dots represent the infrared sources and circles 
mark IRAC sources with X-ray counterparts. We overplot the 
classification of objects dominated by envelope emission (trian- 
gles) and stars dominated by disk emission (crosses). 

of the nebula, according to the criteria discussed in Sect. I4.1.1I 
we therefore consider the remaining 227 sources as candidate 
members of the cluster. 

In Fig(6j we report the color-color diagram of objects in 
the ACIS field of view (without infrared candidate extragalactic 
sources): among the 227 candidate members with X-ray emis- 
sion, we distinguish 116 stars, whose infrared emission is domi- 
nated by disk, 1 object with a strong envelope emission and 110 
sources in the color-color region corresponding to photospheres 
of normal stars. We consider these 110 sources to be Class III 
members of the cluster. 

4.2. Additional YSOs identified via other diagrams 

In the previous sections, we limited our analysis to the stars 
detected in all four IRAC bands, to find the contaminants 
and classify the pre-main sequence stars in our field of view. 
Nevertheless, the above selection criterion strongly limits the 
number of objects, mainly because of the lower sensitivity of 
the 5.8 and 8.0 fim channels: we found in these channels ~ 10% 
of the detections found in the 3.6 or 4.5 fim channels (cf Table 
13. 

Considering the 6768 stars detected in the 3.6 and 4.5 fim 
channels that are not detected in all the four IRAC bands, we 
found 704 stars that have color [3.6] - [4.5] > 0.3 and sharp 
value compatible with that of a point-like source. In the ACIS 
field of view, where there are 3890 infrared detections and 408 
stars with X-ray emission, we found 459 stars with [3.6] - [4.5] 
color excesses, 49 of which present X-ray emission. Because of 
their color typical of PMS stars with disks, these 459 stars are 
candidate members of the cluster, but we are unable to distin- 
guish between protostars and Classical T-Tauri stars because of 
the lack of IRAC color-color information. In Fig. [7] we plot 3.6 



1.0 -0.5 0.0 0.5 1.0 1.5 2.0 

[3.6H4.5] 

Fig. 7. Color-magnitude diagram for stars detected in the 3.6 and 
4.5 fim bands but not in all four IRAC bands. The squares rep- 
resent objects with X-ray emission. The 459 stars with [3.6] - 
[4.5] > 0.3 and sharp parameter compatible with a point-like 
object are candidate members of NGC1893 with infrared ex- 
cesses. 

fim magnitude versus [3.6] - [4.5] color for stars detected in the 
ACIS field of view that have emission in the 3.6 and 4.5 fim 
channels but were not detected in all four IRAC channels, and 
overplot stars with X-ray emission. The main problem related 
to these objects is that we are unable to distinguish them from 
extragalactic contaminants. For this aim, optical spectroscopy 
could be useful to find young objects by means of measures of 
lithium abundance and H a and Ca II lines (or, more basically, 
to distinguish between young stars at zero redshift, and AGN at 
redshifts up to a fe w tenths). 

Other authors (iGutermuth et all l2007t iLucas et all I2007I) 
use also the [[K] - [3.6]]o versus [[3.6] - [4.5]]o color-color dia- 
gram as a means of secondary identification of sources that lack 
[5.8] or [8.0] fim photometric data; due to the sensitivity limit of 
2MASS and distance of NGC1893, in our sample there are, how- 
ever, only 10 stars of the previous 459 candidate members that 
show emission also in the K band, thus for our purposes these 
alternative methods are not useful for identifying new YSOs. 

In Table |4] we report all the candidate members of 
NGC1893; in the first five columns we report the ID and the 
IRAC magnitudes, while in the sixth and seventh columns we 
report two flags, the first indicating if an objects is a Class 0/1, a 
Class II or a Class III star, the second if a source is in the ACIS 
field and if it is an X-ray source. 

5. Summary and Conclusion 

We have analyzed a set of IRAC observations of NGC1893, a 
young cluster, in the outer part of the Galaxy. 

- By analyzing the IRAC images, we were able to identify in 
the field of view (~ 31.2' x 26.0') 1028 sources emitting in 
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Table 4. Candidate Members of NGC1893. ID is the source 
identification number (as in TableQ), [3.6], [4.5], [5.8], [8.0] are 
the IRAC magnitudes, Class is a flag identifying the class of the 
source (0 for Class 0/1 protostars, 2 for Classical T-Tauri stars 
and 3 for Weak-lined T-Tauri), where applicable, X-ray flag is -1 
for stars out of the ACIS field of view, 1 for X-ray sources, for 
stars without X-ray emission in the ACIS field of view 

. The complete list of sources is available in electronic format. 



ID [3.6] [4.5] [5.8] [8.0] Class X-ray flag 

-1 
-1 
-1 



-1 
-1 




-1 
-1 
-1 





1 



-1 
1 

-1 



all four IRAC bands, 25% of which have not been previously 
detected by 2MASS. 

- By analyzing the color-color diagram, we identified 242 stars 
with disks, and 7 Class 0/1 protostars. Moreover in the ACIS- 
I Chandra field of view, we detected 1 10 Weak lined T Tauri. 
At this stage, the census of the cluster is incomplete: we have 
found 459 objects in the IRAC field that showed an infrared 
excess ([3.6] - [4.5] > 0.3) and were not detected at 5.8 
and 8.0 fim. These sources are candidate members of the 
cluster, but we cannot distinguish them from contaminants. 
Moreover, we found a sample of 472 X-ray sources with no 
infrared counterpart, which is probably dominated by extra- 
galactic sources. 

- Although the catalogue of members is incomplete, we have 
detected 359 members, indicating that NGC1893 is quite 
rich, with intense star forming activity despite the "unfavor- 
able" environmental conditions in the outer Galaxy. 

The near and mid-infrared data are unsuitable for estimating 
age and masses of the stars, because the presence of disks af- 
fects the magnitudes. To place our survey in context with other 
IRAC surveys of star forming regions, it is, however, useful to 
derive a rough estimate of the mass range for our census of 
NGC1893. To achieve this, we assume that the J magnitude 
provided by 2MASS is entirely photospheric, for our Class III 
stars. Our brightest Class III star has J ~ 8.28. In order to 
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compare the J luminosity with models, we need to obtain the 
absolute magnitude value that depends on the supposed dis- 
tance ^nd_absorptioji 1 Usmg_d2stence and ab sorpt ion measured 
fromlTapia et al] d 19911) . iMarco et alJ d200lh . and ISharma et al] 
(I2007ir we converted the J magnitudes in absolute magnitudes. 
C omparing our calcula ted values with that obtained from models 
bv lMarigo et alJ d2008l) . we found that the highest mass is in the 
range ~ 28 - 46M . To estimate the lower mass limit, we con- 
sidered that our faintest Class III has [3.6] = 14.6, and converted 
it to absolute magnitude in the same way as for the J magnitude. 
Comparing our va lue with the L magnitudes of PMS models by 
ISiessetalJ d200dt) relative to the ages calculated by previous au- 
thors (1-4 Myr), we estimated that the lowest mass weak-lined 
T-Tauri star is in the range ~ 1.4 - 3M Q . 

We found that the 3% of our stars with IR excesses are Class 
0/1 protos tars. Thi s value is similar to the 2% value found in 
IC348 dLada et al.L 120061) . implying that our cluster does not 
show any peculiarity regarding the ongoing star formation, if 
compared with clusters of the same age in the inner Galaxy. 

The fraction of stars with primordial disks that we measure 
in NGC1893 is about 67%. This value is an upper limit because 
our selection of diskless stars is based on the IR observations, 
therefore, we are able to detect stars with IR excesses to lower 
mass limits than for stars without disks. Considering that clusters 
of similar ages (~ 3 — 4 Myr) may show fractions of disks in 
the range ~ 50-61% wit h significant uncertainty dLada et all 
120061: iLuhman et al.Ll2006h . our upper limit compares well with 
disk fractions measured in the nearby SFRs. 
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